Superconducting vortex profile from fixed point measurements 
The Lazy Fisherman tunnelling microscopy method 
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We introduce a mode of operation for studying the vortex phase in superconductors using scanning 
tunnelling microscopy (STM). While in the conventional STM method, the tip is scanned over a 
sample in which a fixed vortex pattern is prepared, in our "Lazy Fisherman" method the STM tip is 
kept fixed at a selected location while the vortices are being moved by varying the applied magnetic 
field. By continuously acquiring the local tunnelling conductance spectra, dI/dV(V), we detect the 
changes in the local density of states under the tip due to the vortex motion. With no need for 
scanning, the method permits one to extend the study of vortices to samples in which scanning is 
difficult or even impossible due to surface nonuniformity and allows one to study vortex dynamics. 
Using a statistical analysis of the spectra, we reconstruct the single vortex zero bias conductance 
profile. We apply the method to the c-axis face of an MgB2 single crystal sample and obtain a 
vortex profile with a coherence length, £ of 57±2 nm. 

PACS numbers: 74.50.+r 74.70.Ad 74.25.Ha 



In type II superconductors the magnetic field pene- 
trates the superconducting (SC) material, above a cer- 
tain field (Hci), in the form of quantized flux bundles, 
known as vortices, each vortex carrying a flux quantum 
of o =hc/e Q. The vortices arrange themselves in a lat- 
tice and the SC order parameter (OP), A(r), becomes 
spatially inhomogeneous, having a zero value in the cen- 
ter of each vortex and rising to its maximum value in be- 
tween the vortices. Several experimental methods were 
developed for studying the vortex lattice. Most meth- 
ods exploit the spatial variations in the value of the lo- 
cal magnetic field and are therefore sensitive to changes 
occurring over a length scale set by the magnetic field 
jenetration length, A . These include: Bitter decoration 
2j , magneto-optical microscopy Q , local hall microscopy 
4(, electron interference microscopy [{| and squid mi- 
croscopy 0. A different approach was developed using 
a scanning tunnelling microscope (STM)0. The STM 
is used to map the electronic density of states (DOS) as 
a function of position by scanning its tip over the sam- 
ple. This enables one to view the vortex lattice based on 
the observed changes in the value of the Zero Bias Con- 
ductance (ZBC) which manifest the changes in value of 
A(r). Therefore the method is sensitive to changes oc- 
curring on a length scale set by the SC coherence length, 
£. This, together with the STM's ability to obtain atomic 
resolution in the topographic mode, results in images of 
extremely high spatial resolution. The method had been 
successfully used in a number of materials such as NbSe2 
0, BSCCO @, YBCO0, LuNi 2 B 2 c[I3 and MgB 2 [I2. 
A main drawback of this method is that its use is limited 
to highly flat samples. Moreover the relatively large time 
needed for obtaining a single image (several minutes for 
a picture obtained at a single voltage up to hours for the 




FIG. 1: a) A schematic drawing of the experimental method 
showing the STM tip (conic shape) fixed at a chosen position 
and the vortex motion (marked by arrows) . b) A color coded 
3D plot showing 30 normalized conductance spectra as a func- 
tion of voltage and time exhibiting the variations in the local 
DOS under the tip induced by vortex motion. The magnetic 
field time dependence inducing the vortex motion is given by 
the thick line 



full spectra) limits its use in the study of vortex dynam- 
ics. 

In this letter we show how to extend the use of STM to 
samples in which scanning is difficult or even impossible 
by using a different mode of operation which we name 
the lazy fisherman method (LFM) [b| . In LFM the dif- 



2 



ficulties involved in scanning are avoided by keeping the 
STM tip at a fixed position in which at zero field a clear 
SC spectra is observed. The different points within the 
vortex lattice are accessed by creating a vortex motion, 
induced by slowly changing the value of the applied field, 
as illustrated in Fig. la. We show how by statistically an- 
alyzing the data measured at a single fixed point one 
can reconstruct the vortex profile and thus extract the 
SC coherence length. We apply the LFM to a c-axis ori- 
ented surface of a single crystal MgB 2 sample. We find 
that the measured data represents well a statistical en- 
semble sampled uniformly over the vortex unit cell. We 
reconstruct the vortex profile in MgB2 and find the SC 
coherence length to be £=57±2 nm in close agreement 
with the value of 51 nm obtained by Eskildsen et al.[l2j 
using the conventional STM scanning method. 

Single crystals of MgB 2 were grown by a high pres- 
sure method described in Ref. 14] . The experiments 
were carried out on a low temperature Omicron STM. 
The tunneling junctions were achieved by approaching 
mechanically cut Pt/Ir tips to the c-axis oriented sur- 
face of the crystal. Due to the peculiar electronic struc- 
ture of the material, c-axis tunnelling allows to probe 
mainly the quasiparticle density of states in the 7r-band 
|l2l llq . The superconducting critical temperature was 
determined locally by measuring the evolution of the tun- 
neling conductance spectra as a function of temperature, 
and was found to be 38.5 K. The field dependence of the 
tunneling conductance spectra dI(V)/dV was obtained 
by fixing the STM tip in a selected location and continu- 
ously measuring local I—V tunnelling curves while slowly 
sweeping the magnetic field, the field being parallel to the 
c-axis. At low temperature (T = 6.6 K), the field was in- 
creased from zero up to the maximum value of 2.3 T at 
a rate of 0.15 T/min. The acquisition time for a single 
spectrum (-20 mV, 20 mV) was 60 msec. However the 
voltage range in which the DOS is effected by the vortex, 
±A ss ± 2 mV), is scanned in only 6 msec. During such 
a short time the field may be considered constant for an 
individual spectra, as AB ~ 10~ 5 T, thus AB/B « 1. 
As our measured spectra are symmetrical with respect to 
zero bias we conclude that, during the time required to 
obtain a single spectrum, the vortex motion is also negli- 
gible. However subsequently measured spectra do differ 
as can be seen in Fig. lb which shows an example of 30 
consecutively acquired spectra obtained in the field range 
of 0.244 T to 0.250 T. One clearly sees the non mono- 
tonic changes in the local DOS under the tip, due to the 
vortex motion. The spectra range from a clearly gapped 
form showing shoulders at the gap edges (reflecting the 
DOS far from the vortex core) to a completely flat form 
(reflecting the DOS at the vortex core). 

A commonly used parameter in conventional STM vor- 
tex mapping is the value of the normalized ZBC, a = 
[dI/dV(V = 0)}/[dI/dV(V » A)]. In the normal state 
this value is equal to 1 while in the SC state it is strongly 
reduced due to an opening of a gap in the DOS around 
the Fermi level. Thus, one obtains the vortex image by 
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FIG. 2: Field evolution of the Normalized ZBC histograms: 
T, 0.24 T, 0.34 T, 0.41 T, 0.49 T and 2.08 T. Histograms 
obtained at intermediate field values (0.24 T - 0.49 T) show 
an asymmetric distribution peaked at or slightly above the 
minimum, reflecting the shape of the vortex profile. Inset 
shows the ZBC values obtained for 0T <H<0.5 T. Rectangles 
mark the ranges used for obtaining the histograms shown in 
the main figure 



mapping the value of the ZBC as a function of position. 
The ZBC obtains its maximum value at the vortex core, 
gradually decreases as a function of the distance from 
the vortex core and reaches its minimum value at the 
point of equidistance from the vortices (see for example 
[12 )■ Measuring the value of the ZBC as a function of 
the distance from the vortex core allows one to obtain 
the exact vortex profile and thus to estimate the SC co- 
herence length, £. 

As in our method one can not directly determine the 
distance between the vortex core and the STM's tip we 
use a statistical method in order to study the distribution 
of the ZBC values and to reconstruct the vortex profile. 
As a first step we have created histograms of the ZBC 
values acquired over 10 sec intervals. The time interval 
was chosen to be on one hand short enough to assure that 
the inter-vortex distance does not significantly change 
(Ad/d < 6%) and on the other hand to allow enough 
data points for statistical analysis. 

The inset of Fig. 2 shows the ZBC values as a func- 
tion of field for 0<H<0.5 T. The field in which the vor- 
tex start to pass under the tip is clearly seen as a sharp 
jump in the values of the ZBC (see arrow in the inset). 
As this value is always larger than Hci due to surface 
barrier effects we conclude that Hci|| c < 0.12 ± 0.03T 
in agreement with the value of 0.1 T found by Lyard 
et al.|l6j by magnetization measurements. The resulting 
histograms at fields of T, 0.24 T, 0.34 T, 0.41 T, 0.49 
T and 2.08 T are shown in Fig. 2. Each histogram rep- 
resents roughly 150 ZBC values obtained in a field range 
of 0.03 T. The zero field histogram is narrow, symmet- 
ric and as expected shows no sign of spread in the ZBC 
values due to vortices. At higher fields (H>0.1 T) the his- 
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FIG. 3: vortex profile, c(r) (circles) with a theoretical de- 
pendence given by l-tanh(r/£)(line). Best fit is given by 
£=57±2nm. Inset shows the Renormalized ZBC values his- 
togram used to calculate the vortex profile 

tograms become highly asymmetric. The probability is 
peaked near the minimum conductance value and decays 
considerably towards the maximum value. Looking at 
the set of histograms one can see an increase of the mini- 
mum value as the field increases, reflecting the increasing 
overlap between the vortices. In the case of MgB2 the 
exact field dependence of this increase is rather compli- 
cated due to coupling between the two bands. We will 
discuss the properties related to the unique two band na- 
ture of MgB2 in a separate publication [L7J. Finally at 
high fields (H>2 T) the histogram becomes symmetric 
again with the modulation due to the vortices being so 
small that it is almost entirely washed out by the exper- 
imental noise. As no gap can be detected in any of our 
spectra for H>2.3T we deduce H c2 || c (T = 6.6K) > 2.3 
a value comparable to that measured by Lyard et al.|l6| 
of 2.6 T. 

Using the ZBC histograms we are able to re- 
construct, with no free parameters, the vortex pro- 
file, ZBC(r)(=<j(r)), where r is the distance from 
the vortex core. Fig. 3 shows the results calcu- 
lated from the experimental histogram obtained in 
the field range of 0.24-0.27 T (see inset of Fig.3). 
The abscissa which is the distance from the vortex 
core is calculated directly from the histogram using 

r(fJ n )/r(a min )=^Jj2i= n +i-,N p { a %^ a i + Sa )> where 8a 
is the bin size. This expression is exact for a vortex 
profile of cylindrical symmetry. In such a case the proba- 



bility of obtaining a ZBC value between o\ and a 2 which 
is general given by P(ai,a 2 )oc J d9(r 2 (a2,0) — r 2 (cr 1 ,6')) 
reduces to P(<7i,(T2)oc (r 2 (o2) — r 2 (<7i) . This approxima- 
tion is valid as long as the inter-vortex distance is large 
in comparison to the coherence length (£ <d) and that 
the material is in the dirty limit (£ > £), which is the case 
111 MgB 2 [T2. The ZBC is presented on a renormalized 
scale (cr(ff) - a{H = 0))/(l - a(H = 0)). This allows 
one to separate the contribution of DOS changes to the 
ZBC from that of finite temperature effects and inelastic 
tunnelling process. The solid line in the figure is a best fit 
curve to the data obtained using c(r)=l — tanh(r/£). The 
best fit is obtained for a value of ^/r(cr m j n )=0.99±0.03. 
As the minimum ZBC value is obtained (for a triangu- 
lar lattice) at r(cr mm ) = d/ v / 3=59 nm (for H = 0.24 T), 
where d is the inter vortex distance given by d ~ 50 
nm/VH, this results in £ = 57 ±2 nm. 

To conclude we have presented a method which allows 
to extend the use of STM for vortex study to cases where 
scanning is difficult or impossible and to cases showing 
fast vortex dynamics. Using the method one can measure 
several basic parameters describing the SC vortex state 
such as the field of first vortex penetration, the SC co- 
herence length £ and the upper critical field HC2 . At the 
same time our method maintains the scanning capabili- 
ties of the STM. Thus, it allows repeating the measure- 
ment at different fixed points in the sample and relating 
the observed changes in vortex behavior to the position 
in the sample. This could be of interest in the vicinity of 
impurities, pinning centers and grain boundaries. Addi- 
tionally the LFM is suitable in all type of samples (rough 
and flat) for studying faster vortex dynamics (in com- 
parison to conventional scanning) being limited only by 
the time needed to measure a single spectrum, ~lms. A 
variant of the method can be used for studying faster vor- 
tex dynamics by measuring the tunnelling current (with 
open feedback loop) at a voltage of eV < A. As vortex 
motion is the source of dissipation in the SC state, LFM 
being a local fast tool to study the SC parameters re- 
lated to the vortex lattice and the vortex dynamics could 
of great use in the development of SC applications. 
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